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Abstract 

The large value of #13 recently measured by rector and accelerator experiments opens un- 
precedented opportunities for precision oscillation physics. In this paper, we reconsider the 
physics reach of medium baseline superbeams. For #13 ~ 9° we show that facilities at medium 
baselines - i.e. L ~ O(1000km) - remain optimal for the study of CP violation in the leptonic 
sector, although their ultimate precision strongly depends on experimental systematics. This 
is demonstrated in particular for facilities of practical interest in Europe: a CERN to Gran 
Sasso and CERN to Phyasalmi beam based on the present SPS and on new high power 
50 GeV proton driver. Due to the large value of #13 , spectral information can be employed 
at medium baselines to resolve the sign ambiguity and determine the neutrino mass hierarchy. 
However, longer baselines, where matter effects dominate the — > v e transition, can achieve 
much stronger sensitivity to sign(Am 2 ) even at moderate exposures. 

1 Introduction 

After an experimental search lasted more than a decade [2 El El 2] , we know quite precisely 
the size of the mixing angle between the first and third neutrino family (#13) El El |S]. 
Such value turned out to be extremely large - in fact, very close to previous limits set by 
CHOOZ [9] and Palo Verde [10] . A large value of #13 not only constraints models for neutrino 
mass and mixing but, even more, opens up the possibility to perform precision physics in 
the leptonic sector with artificial neutrino beams in a way that resembles what was done for 
quark mixing in the B-factory era [12] . In the past, several experimental strategies have been 
considered |1 31 114) . corresponding to possible values of 613 . The most extreme, addressing 
sin 2 2^13 as small as 10 -4 , required novel acceleration techniques to achieve unprecedented 
neutrino beam intensities and purity. The results from T2K, Daya-Bay, RENO and Double- 
Chooz allow to reconsider these scenarios, seeking for less challenging experimental setups to 
establish CP violation in the leptonic sector and determine the neutrino mass pattern ("mass 
hierarchy"). The experimental proposals will take advantage of the large i/ M — > oscillation 
probability due to #13 = 8.9° ± 0.4° [15) . likely relieving the constraints on the detector mass 
and accelerator power. In this paper we address some of these issues, firstly on general ground 
and, then, driven by practical considerations. On general ground we discuss to what extent 
the constraints on beam power, purity and systematic uncertainty can limit our capability to 
establish CP violation in the leptonic sector using Superbeams [161 1171 1181 1191 1201 I21| . We 
consider the optimal baseline to search for CP violation at the atmospheric scale, i.e. a medium 
baseline (L ~ 1000 km) with neutrino energies of E ~ 1 GeV and a far detector capable of 
reconstructing quasi-elastic, deep-inelastic and resonance final states (liquid argon TPC). It is 
well known that such baseline represents a poor choice for the simultaneous determination of 
CP violation (CPV) and mass hierarchy, which could be better achieved with longer baselines 
and wide- band neutrino beams [221 123] . Hence, in the context of medium baselines we evaluate 
the deterioration of the sensitivity on CPV due to the ignorance on the neutrino mass hierarchy. 
Similarly, we compare the CP reach of the facility when the far detector is located at a baseline 



optimal for the determination of mass hierarchy through the exploitation of matter effects 
[L 3> 1000 km). The results of such general study are relevant in setting experimental issues 
and to establish a European strategy for precision measurements in the large 6*13 era. 

First, we discuss to what extent the Gran Sasso underground Laboratories could be consid- 
ered as a viable candidate to host the facility for CP violation. As an alternative option [24] . 
we detail the physics reach of a new shallow depth laboratory able to host significantly larger 
detector masses than what can be accommodated in the present experimental halls of LNGS. 
Finally, we compare the CP reach of these facilities with a multipurpose detector [251126] located 
at longer baselines, i.e. from ~1000 up to 2290 km (CERN-to-Phyasalmi [27] [28l [29] [30]), in 
order to establish to what extent the choice of a baseline better suited for the mass hierarchy 
can affect the CP reach. All comparisons are done as a function of the intensity of the neu- 
trino source, detector mass and systematic uncertainties. When needed, beam parameters are 
re-optimized ab initio (i.e. at the target-horn level) to achieve best performance for a given 
configuration and allow a fair comparison among various options. 



2 Facilities at medium baselines 

All current long-baseline experiments are in the medium baseline range (~ 1000 km) either 
because they are tuned to be at the peak of the oscillation probability at the atmospheric 
scale (MINOS: 730 km, T2K: 293 km, NOVA: 810 km) or because they maximize the event 
rate for f M — > v T transitions (OPERA: 730 km). At these baselines, matter density can be 
safely considered constant and — > v e transition probabilities are commonly expressed by a 
perturbative expansion on sin26*i3 and a = Am^i/Am^. In the following, P{y^ — > u e ) is 
approximated as |31j : 



P(^,, -> u e ) ~ sin 2 20i3 sin 2 #23 
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In this formula A = Am^ 1 L/(4:E) and the terms contributing to the Jarlskog invariant are split 
into the small parameter sin 26*13, the 0(1) term £ = cos 6*13 sin26*i2 sin26*23 and the CP term 
sin <5; A = 2\/2Gpn e E / Am^i with Gf the Fermi coupling constant and n e the electron density 
in matter. Note that the sign of A depends on the sign of Am 2 i which is positive (negative) for 
normal (inverted) hierarchy of neutrino masses. 

The large size of 6*13 changes the phenomenology of the medium baseline. In the past it 
was commonly believed sin 26*13 to be at most of the size of a, i.e. a large suppression of 
the first term in Eq[l]was expected. We know now that sin26*i3 overwhelms a by one order of 
magnitude (~0.3 versus ~0.03), which makes the determination of the CP-blind term 0\ crucial 
for precision oscillation physics. Since Ox is large and depends on the mass hierarchy through 
A, the ignorance on the sign of Am^i plays a role even at medium baselines. 

For sake of illustration, Fig. [I] shows the oscillation probability V(y^ — \ v e ) versus neutrino 
energy for different choices of the mass hierarchy and the 8c p phase for a baseline of 730 km and 
sin 2 26*13 = 0.092. It is evident that effects ruled by the mass hierarchy (matter effects) and the 
CP phase are of comparable size. At first order they introduce a change in the normalization so 
that the control of systematic errors becomes crucial. The second oscillation maximum, that can 
be used to disentangle CP and matter effects, appears at an energy of about 500 MeV. The full 
exploitation of spectral information down to a few hundreds of MeV points toward a detector 



with high granularity, energy resolution and efficiency for low-energy electron neutrinos. That 
justifies the consensus on the use of liquid argon (LAr) detectors as far detectors for precision 
experiments at medium baselines (see Sec.[3]|. 

In the following, three configurations will be detailed. All of them exploit a neutrino beam 
from CERN to the Gran Sasso area and are based on a 

1. on-axis detector inside the existing underground laboratory (up to 10 kt mass) (ON A) 

2. 7 km off-axis detector (up to 100 kt mass) located in a new shallow depth site (OA7) 

3. as above but at 10 km off-axis (OA10) 

In particular, OA7 has been studied in the past [25 for the measurement of #13 although the 
mean neutrino energy is quite far from the oscillation peak; on the other hand, the spectrum of 
OA10 better matches the region of interest for CP violation at the price of reduced statistics. 
The first option (ONA) is clearly of great practical interest since it is completely based on 
existing underground facilities. However, the limited size of the LNGS halls strongly constraints 
the maximum detector weight [32]. The ICARUS T600 module [331 134) . which is currently in 
data taking in the Hall B of LNGS (110 m in length), has a dimension of 3.9 x 4.3 x 19.6 
m 3 , corresponding to a liquid argon mass of 0.735 kt (0.476 kt fiducial). The ICARUS T1200 
module 35 proposed in 2001 for the CNGS was based on a tank of 10.3 x 10.3 x 21 m 3 , totaling 
1.47 kt (0.952 kt fiducial). Assuming an array of four T1200 modules located in Hall B would 
allow the installation of a 4 kt fiducial mass detector at Gran Sasso. In addition, a novel design 
recently considered by the ICARUS Collaboration can increase the total available mass in Hall 
B up to 7.5 kt [5S]. In conclusion, even assuming the availability of extra space at LNGS, the 
Gran Sasso laboratories could not host a detector with mass larger than 10 kton. In this study 
10 kton is considered the largest explotable mass for the ONA option. 

Similarly, two classes of options have been considered for the neutrino source (see Sec. 131. 
The former, which is based on the existing CERN-SPS, employs at most facilities that are 
already available; the second one, which considers a new CERN-based proton driver at 50 GeV, 
allows for higher flexibility in the choice of the beam configuration and better performance on 
a longer timescale. 
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Figure 1: Viy^ — > v e ) (left) and Viy^ — > v e ) (right) vs neutrino energy for different choices of the 
mass hierarchy and 5 for a 730 km baseline and sin 2 2# 13 = 0.092. 




3 Beamline and detector simulation 

The beam parameter optimization for a novel 50 GeV proton driver is based on the results 
obtained in [37] in the framework of the LAGUNA design study. The most relevant parameters 
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Figure 2: Focusing optics for the 50 (top) and 400 GeV (bottom) proton energy beamlinc. 



are the distances among the target, the upstream focusing horn and the downstream focusing 
horn ("reflector"), the target length and the length and width of the decay tunnel. Following 
the approach of Ref. [37], the parameters, together with the geometry of the horn and reflectors 
were optimized using the sensitivity to #13 as a figure of merit, instead of resorting to inter- 
mediate observables as e.g. the number of unoscillated neutrinos. This optimization is also 
appropriate for the CP reach, especially if the detector technology allows the exploitation of 
spectral information. The optimal configuration for a ~ 700 km baseline (CERN-to-Gran Sasso) 
is shown in the top drawing of Fig. [2] The decay tunnel is 90 m long and the optimal radius is 
2.2 m. The target consists of a 1 m long graphite rod. Primary interactions in the target were 
simulated with the GEANT4 [38] QGSP hadronic package and secondaries were traced down to 
the decay tunnel. The #13 sensitivity was computed using the GLoBES 3.1.11 package [391 14U] 
and assuming 100 kton LAr far detector (see below). 

The neutrino beam from the existing SPS required a dedicated optimization since the current 
CNGS configuration (on-axis, (E v ) ~ 17 GeV) is not appropriate for CP violation studies. An 
approach similar to the optimization of the 50 GeV driver was followed. A pre-scan of the beam 
parameters was carried out using a fast simulation based on the BMPT parameterization [41] , 
This approach is particularly rewarding since the BMPT formulas are mostly based on data col- 
lected with 400 GeV/c and 450 GeV/c protons at CERN-SPS. However, fluxes in the proximity 
of the optimal parameters were fully simulated with GEANT4. The optimised configuration is 
shown in the bottom drawing of Fig. [2] The tunnel geometry is the one of CNGS with a length 
of 1000 m and a radius of 1.225 m. A 1 m long graphite target was used. The corresponding 
off-axis event rates in the optimal configuration are shown in Fig. [3] together with the rates for 
the ONA configuration ("50 GeV ON"). All rates are normalized to the proton energy. 

The advantage of a dedicated proton driver at 50 GeV is quite evident since the region where 
most of the CP related information reside is located below 2.5 GeV. Moreover, at these base- 
lines the second oscillation maximum is inaccessible for any configuration based on a 400 GeV 
accelerator. Similar considerations hold for a low energy (10 GeV optics) re-optimization of 
CNGS ("CNGS LE on-ax" in Fig. [3), as the one considered in Ref. [42]. 

Liquid Argon detectors are the most promising technology to address precision oscillation 
physics. The main advantages of a LAr TPC in v e appearance are very high efficiency both for 
quasi-elastic (80%) and deep-inelastic (90%) interactions combined with a superior NC rejection 
power. In the following, we considered a contamination of NC due to 7r — > e misidentification 
not exceeding 0.1% of the CC rate. In the occurrence of QE interaction, the neutrino energy 
can be fully reconstructed by the lepton energy and direction, since the direction of the incoming 
neutrino is known in advance. LAr detectors [33] [34J |33] are capable to reconstruct E v with a 
resolution mostly dominated by the electron energy resolution: oe v /E v ~ 0.05/V 'E v , E v being 
expressed in GeV. On the other hand, energy resolution for DlS-i^ interaction is driven by the 
resolution on the hadronic system. In LAr, the latter amounts to aE h /E h ~ 0.2/ ^/ E h (GeV). 
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Figure 3: v^ c event rates for the two off-axis (OA) options considered in the text ("400 GeV OA 
7-km" and "400 GeV OA 10-km") and for a dedicated on-axis 50 GeV proton driver ("50 GeV ON"). 
"CNGS LE on-ax" refers to the low-energy (10 GeV optics) re-optimization of CNGS considered in 
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In this study, the energy resolution and efficiency for v £ and i> e were implemented smearing the 
final state momenta of the electrons and hadrons. Interactions were simulated using the GENIE 
Monte Carlo generator |44) and the corresponding migration matrices were implemented in the 
detector description of GLoBES. The smearing matrices were calculated for and z/ e separately. 
They are shown for v e events in Fig. [4] (left). In the simulation, a 50 MeV bin width was used. 
Similarly, the efficiency as a function of energy is shown in Fig. [4] (right). 

4 Results 

As mentioned above, three possible options were considered in this study: an on-axis detector 
with limited fiducial mass (ONA), a 7 km (OA7) and 10 km (OA10) off-axis setup leveraging the 
existing 400 GeV accelerator. In Fig. [HJthe v^ c appearance spectra convoluted with detector 
effects for the ONA, OA7 and OA10 options are shown. For the on-axis option, a 50 GeV proton 
driver with a 10 kt detector was assumed. The power of the proton driver corresponds to 3 • 
10 21 pot/y (2.4 MW). For the off-axis options, a 20 kton detector and 1.2 ■ 10 20 pot/y (0.77 
MW) was considered. Clearly, for the 400 GeV proton driver option the power of the driver is 
constrained by the limitations of the SPS-based neutrino beam, while for the 50 GeV facility, a 
dedicated multi-MW machine can be envisioned. 

The performance of the different options were compared defining as figure of merit the CP 
violation discovery potential at 3a. Following [13U45] , we define the "CP coverage" at 3<j as the 
fraction of possible (true) values of the CP phase 8 where the CP conserving hypothesis (8 = 0, n) 
has a p-value smaller than 0.0027. It corresponds to x-min > 9 for 8 = 0,7T (1 d.o.f.). In each 
case, we considered separately the possibility that the mass hierarchy is known by the time the 
facility starts data taking. If the mass hierarchy is not known, the p-value of the null hypothesis 
is computed as the Xmm of the following combinations: (8 — 0, Am^ > 0),(<5 = ir, Am| 2 > 0), 
{8 = 0, Am| 2 < 0), {8 = 7T, Am| 2 < 0). 

The 3a CP coverage in the on-axis facility based on the 50 GeV proton driver is shown in 
Fig. [6] for 5 years of running with v and 5 years with v. The corresponding results are shown in 
Figs!]7]and [5] for the 7 and 10 km off- axis configurations based on the 400 GeV proton driver. In 
each plot we consider normal and inverted hierarchy, assuming this information to be available 
or not (color codes) at time of running. The assumed systematic error on flux normalization 





Figure 5: v e appearance signal at 730 km for 5 years of data taking, sin 2 26> 13 = 0.092, and normal 
hierarchy. Upper Left: ONA option, 10 kt <g> 3 • 10 21 pot/y (2.4 MW). Upper Right: OA10 option, 
20 kt ® 1.2 ■ 10 20 pot/y (0.77 MW). Lower Left: OA7 option , 20 kt ® 1.2 • 10 20 pot/y (0.77 MW). 
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Figure 6: CP coverage at 3ct for 5 years of v and 5 years of v running with the on-axis configuration, 
a 50 GeV proton driver and a systematic error on flux normalization of 5%. We consider normal 
and inverted hierarchy assuming this information to be available or not (color codes). 



amounts to 5 %. 

A comparison of the three options in the case of normal hierarchy (assumed to be unknown) 
is plotted in Fig. [9] 

Some interesting considerations hold either when the results are considered on general ground 
for any medium baseline facility or when specific CERN-based facilities are envisioned. Firstly, 
it is clear that the limitation due to the sign degeneracy |46i 147| , i.e. the deterioration of the 
CPV coverage due to the ignorance of the mass hierarchy plays a major role only for limited 
exposures or detector mass and it is not an intrinsic limitation of medium baselines for such 
high values of #13. Large exposures, however, are mandatory to exploit spectral information 
and thus to remove the sign ambiguity at L ~ 730 km. 

The constraints on the minimum exposure is particularly severe when we consider realistic 
CERN-based facilities. In the vertical lines of Figs. [6] we indicated the CP coverages of a facility 
whose far detector is hosted in the underground halls of LNGS (maximum LAr mass: 10 kton). 
It is apparent that LNGS is not appropriate as a far detector site unless a dedicated proton 
driver well exceeding 2 MW power can be built at CERN. Similar considerations hold for the 
off-axis facilities that leverage the infrastructure on which the CNGS is built (Fig 7|8 1. In this 
case, a new shallow laboratory can overcome the limitation on the detector mass but the low 
power of the driver constraints the CP coverage well below 50%. Again, the limited statistics 
reduce the resolution power of spectral information to lift the sign ambiguity, so that prior 
knowledge of the mass hierarchy increases the CP coverage up to a factor 2. 

At large exposures the coverage quickly saturates at ~70% (see Fig[9|. This is an intrinsic 
feature of Superbeams for large #13 where the CP reach is dominated by the systematics on the 
flux normalization. That is shown in Fig. |10| where systematics are varied from 1% to 10%. In 
this range, the coverage at saturation decreases from 86% to 60%. This result confirms the key 
role played by the near detector and by the knowledge of the cross sections and detector efficiency 
when performing precision physics with Superbeams [48]. In fact, a control of the systematics 
at the few-percent level remains very challenging, being still at the level of 10.3% for T2K. 
Clearly, this issue makes the construction of a near detector mandatory for CPV searches. This 
is particularly relevant for CERN-based facilities, where the installation of a near detector poses 
additional difficulties from the point of view of engineering and civil infrastructure. 

Finally, it is interesting to compare the CP reach of medium baselines with facilities that 
exploit longer distances to establish the mass hierarchy through the observation of matter effects. 




Figure 7: CP coverage at 3<r for 5 years of v and 5 years of v running with the off-axis 7 km 
configuration, a 400 GeV proton driver and a systematic error on flux normalization of 5%. We 
consider normal and inverted hierarchy assuming this information to be available or not (color 
codes) . 




Figure 8: CP coverage at 3er for 5 years of v and 5 years of v running with the off-axis 10 km 
configuration, a 400 GeV proton driver and a systematic error on flux normalization of 5%. We 
consider normal and inverted hierarchy assuming this information to be available or not (color 
codes). 
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Figure 9: Comparison of the CP coverage at 3<r for 5 years of v and 5 years of v running with the 
on-axis 50 GeV configuration, off-axis 7 and 10 km configuration and 400 GeV proton driver, and 
a systematic error on flux normalization of 5 %. Here we only consider normal hierarchy, which is 
assumed to be unknown. 
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Figure 10: Effect of systematics on the CP coverage for the on axis (left) and off-axis 10 km (right) 
configurations. 
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Figure 11: Exposure needed to reach 40% (solid line) and 70% (dashed line) CP coverage as a 
function of the baseline. The red, blue and purple contours are drawn considering the optimal 
beamline (see Sec. [3]) for a 730, 1570 and 2300 km baseline, respectively. 



In Fig.[TT]we show the 70% and 40% CP coverage contour as a function of baseline and exposure 
for three different beamline configurations. The configurations are labeled $730, $1570 and 
$2300: they correspond to beam parameters optimized (see Sec. [3| for L — 730, 1570,2300 km, 
respectively. In addition, Fig. 12 shows the 3<r CP coverage for the CERN-to-Phyasalmi (2290 
km) facility based on a dedicated 50 GeV proton driver. The same facility without the dedicated 
proton driver, i.e. leveraging the existing SPS complex has a CP coverage described in Fig. |13| 
As far as mass hierarchy is known in advance, medium baselines from 700 to 2300 km exhibit 
similar performance to establish CP violation in the leptonic sector. In fact, larger baselines are 
slightly favoured due to the broader neutrino energy range and, hence, due to the increase of 
spectral information. In this framework, the optimal baseline is somehow in between LNGS and 
Phyasalmi, i.e. at L ~ 1800 km (Fig. 11 1. At L ~ 2000 km the CP reach is independent of the 
mass hierarchy because the sign of Am 2 3 can be established during data taking for all values of 
8. This is not the case for facilities at L ~ 730 km, where the mass hierarchy sensitivity hardly 
exceeds the one of the NOVA experiment [49] (see Fig. 141. Here, the mass hierarchy coverage 
is defined as the fraction of possible (true) values of the CP phase 8 where the wrong hierarchy 
hypothesis is disfavoured at 3<r level. 



5 Conclusions 

The remarkable size of #13 makes Superbeams the most straightforward choice for precision 
physics in the neutrino sector. While the determination of the mass hierarchy seems at reach 
for the next generation of facilities, the study of CP violation remains a major experimental 
challenge. Superbeams at medium baselines offer a clean environment to study CPV, provided 
that the mass hierarchy is known by the time of running or the sign ambiguity is resolved 
employing spectral information. However, major technological efforts are required to attain 
appropriate neutrino fluxes. In this paper we reconsidered medium baseline superbeams with 
special emphasis on CERN-based facilities, which have been re-optimized simulating the beam- 
line from target to the decay tunnel. Facilities that are based on the CERN-SPS and exploit 
the CNGS beamline (off-axis 10 km) reach a modest CPV coverage in any realistic configu- 
ration: at most 52% (48%) with a 100 kton detector in a new shallow laboratory, assuming 
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Figure 12: CP coverage for the facilities based on the 50 GeV proton driver (see text for details). 



CP coverage at 3 a (%), 5+5 y, err.sys. = 0.05 

lo 100r ; ; ; ; : 

p _ ^— ^— NH known. OFA-1 Okm 400(3eV, 730 km \ \ \ 

° g0 ~ NH unknown. OFAi-10km 400GeV, 730 km ! ! I 

~ NH known. ONA5pGeV, 22|90 km j j | 




Figure 13: CP coverage for the facilities based on SPS (400 GeV proton driver) 10 km off-axis, and 
on a dedicate driver (50 GeV) at L=2290 km. 
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Figure 14: Mass hierarchy coverage for the 50 GeV on-axis option (left) and the 400 GeV 10 km 
off- axis option. 



the mass hierarchy to be (not) known. Hence, a dedicated high-power driver is mandatory for 
CPV studies. Similar considerations hold at longer baselines. Here we considered specifically 
the CERN-to-Phyasalmi distance, corresponding to 2290 km, and a new on-axis neutrino beam 
from CERN. For a detector mass of 100 kton the CPV coverage is 62% and drops to 28% if the 
LAr mass is lowered to 20 kton. On the other hand, at large baselines, the size of matter effects 
allows for a determination of the mass hierarchy for any value of 8. Hence, the facility does not 
suffer from the sign ambiguity even if the sign of Am^ is unknown a priori. The optimal facility 
to address CPV is an on-axis configuration (see Figs. 9| 10 \ with a dedicated high power proton 
driver. A 100 kton detector hosted in a new underground facility could reach a CPV coverage 
of ~ 70%, fully dominated by systematic uncertainties: 70% coverage for an overall systematic 
uncertainty of 5% and 86% coverage if the systematics can be brought down to the (currently 
unrealistic) value of 1%. Hence, in way similar to what happened with reactor experiments from 
Chooz to Daya Bay, novel techniques must be developed to lower systematics at the percent 
level in order to overcome the O(70%) coverage limit. 
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